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Abstract 

A new non-relativistic expansion in terms of the nucleon's momentum inside nuclear matter 
of the current for isobar electro-excitation from the nucleon is performed. Being exact with 
respect to the transferred energy and momentum, this yields new current operators which re- 
tain important aspects of relativity not taken into account in the traditional non-relativistic 
reductions. The transition current thus obtained differs from the leading order of the traditional 
expansion by simple multiplicative factors. These depend on the momentum and energy transfer 
and can be easily included together with relativistic kinematics in non-relativistic, many-body 
models of isobar electro-excitation in nuclei. The merits of the new current are tested by com- 
paring with the unexpanded electromagnetic nuclear responses in the isobar peak computed in 
a relativistic Fermi gas framework. The sensitivity of the relativistic responses to the isobar's 
magnetic, electric and Coulomb form factors and the finite width of the isobar is analyzed. 
PACS: 25.30.Rw, 14.20.Gk, 24.10.Jv, 24.30.Gd, 13.40.Gp 

Keywords: Nuclear reactions; Inclusive electron scattering; Delta isobar electro-production. Rel- 
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1 Introduction 



The cross section for inclusive electron scattering (e, e') shows a pronounced peak at an 



energy transfer cu ~ \jq 2 + m 2 N — m^, corresponding to the quasi- free interaction with 
the individual nucleons in the nucleus (here ttin = nucleon mass). For high values of the 
momentum transfer q — |q| and higher energy loss it is possible to produce real pions and 
the cross section shows another peak dominated by the resonant production of a A(1232) 



at uj ~ y q 2 + m 2 A — m 2 N , where ira is the A mass [|l| . The width of these peaks is related 
to the Fermi momentum of the nucleons inside the nucleus and, in the case of the A-peak, 
also to the decay width of the A in nuclear matter. Hence for a high enough value of q, 
these two peaks actually overlap and cannot be separated in inclusive experiments |2|. |3[ . 
Thus the response in the region above the quasi-elastic peak contains information about 
the nucleon's excited states and their change due to the nuclear medium. 

Since the electro-excitation of the A requires high energy and momentum transfers, a 
relativistic treatment of the reaction is needed. Recently, several many-body calculations, 
both in nuclear matter and finite nuclei, have been performed in this region @, ||, ||; all 
of these calculations are non-relativistic in nature, although some relativistic corrections 
enter in two of them, including an expansion of the current to order (p/itln) 2 in @] and 
relativistic kinematics in |J. Clearly some of these corrections are inadequate when one 
wishes to go to high momentum transfers, q ~ 1 GeV/c, and there relativistic models 
such as the ones developed in 0, ||, ^ |T(J are more appropriate. However, although these 
last calculations are fully relativistic, they do not include the full N-A vertex. For in- 
stance, in the pioneering calculation by Moniz |7]] the Peccei Lagrangian was used, which 
is only appropriate for computing the transverse response for low momentum transfer 



[fll| , |12| , pT3| . In other work j8], |9| a more appropriate Ml magnetic transition current was 
used, although the electric E2 and Coulomb C2 excitation amplitudes of the A were not 
included. For years an important program has been pursued to determine more accurately 
the quadrupole C2 and E2 amplitudes in the A region |14[, these being small compared 
with the dominant dipole Ml amplitude. Using polarized photons, the E2/Ml-ratio has 
been measured to be around -3% at resonance |15| . The C2 amplitude however only ap- 
pears in electro-production reactions N(e, e')A. Values of the ratio C2/M1 around —13% 
have been reported in H(e,e'ir°)p experiments at Q 2 = 0.13 (GeV/c) 2 ||16|| . However, 
this value differs with the findings of recent measurements of the transverse-longitudinal 
asymmetry and proton polarization in H(e,e'p)n° reactions |17| . Hence, for the N — > A 
transition our knowledge is still incomplete and has not been possible to undertake a full 
analysis in the sense of the work by Nozawa and Lee [113 of the effect of the C2 and E2 
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form factors in the nuclear A-peak. 

In this paper we perform a new non-relativistic expansion of the electro-excitation 
current of the nucleon, providing an extension of our previous expansion of the electro- 
magnetic nucleon current in powers of r\ = p/m^, in which we retained the full dependence 
on q and u |19[ . Recently the same procedure has been also applied to meson-exchange 
currents (MEC) ||20|| . These currents can be implemented together with relativistic kine- 
matics in standard non-relativistic models of one-particle emission near the quasi-elastic 
peak. In this paper we apply the same procedure to the A electro-excitation current, 
which we develop to leading order in rj, again retaining the full dependence on the energy 
and momentum transfers. We perform this expansion for the magnetic transition current 
Ml, which is the dominant one both in the longitudinal and transverse nuclear responses 
. The resulting current is designed in such a way that it differs from the traditional non- 
relativistic limit simply by having (q, tu)-dependent factors which multiply the traditional 
operators. These corrections, being of leading order in 77, are seen to arise mainly from 
the normalization factor in the Rarita-Schwinger spinor and from the lower-component 
spinology now included in the effective current operator. 

The organization of the work is as follows. In sect. 2 we first develop the analytical 
expressions for the longitudinal and transverse responses in the relativistic Fermi gas 
(RFG) model by using a A-hole approach. We consider the full vertex of Jones and 



Scadron [12| that includes Ml, E2 and C2 A-amplitudes. In deriving expressions for these 
responses we assume a stable A particle and later include its finite width by performing 
a convolution of these responses with a Lorentz distribution. In sect. 3 we perform the 
expansion of the magnetic current to leading order in the momentum of the bound nucleon. 
In sect. 4 we test the validity of the expansion by comparing the exact RFG result with 
a non-relativistic Fermi gas model, using the new current and relativistic kinematics. 
Furthermore, we compare several models of the reaction by using the RFG; in particular, 
we study the differences that arise upon using the Peccei and magnetic Lagrangians, 
and explore the effects of the E2, C2 multipoles and the finite width of the A on the 
longitudinal and transverse response functions. Finally in sect. 5 we draw our main 
conclusions. 



2 General formalism 

We start our discussion by introducing the general formalism, which is based on a rel- 
ativistic treatment of the nuclear currents entering into the calculation of the response 



functions. It is well-known [21] that the longitudinal and transverse response functions 
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can be evaluated as components of the nuclear tensor Wn U , namely 



R L (q,uj) = — 



Q 2 , 



Wt 



oo 



- (w 03 + w 30 ) + -^w 33 

q q A 



W, 



00 



R T (q,uj) 



W u + W, 



22 



(1) 
(2) 



where = (00, q) is the space-like four-momentum carried by the virtual photon and 
the gauge invariance has been exploited in obtaining Eq. dH). We shall compute this 
nuclear tensor in the RFG framework, where nucleons are assumed to move freely inside 
the system with relativistic kinematics, hence being on their mass-shell. 



2.1 Nuclear tensor and response functions for inelastic processes 

We first consider the electro-production of a stable resonance (namely, the A is viewed 
as a particle on its mass-shell) in the RFG - - in this case analytical expressions for 
the response functions are obtained — and later on we shall include corrections due to 
inclusion of the decay width of the A. The RFG nuclear tensor reads (see, for example, 

P) 



W, 



3vrWm 2 A 



N 



Z.3 
hp 



dp 6{k F — p) 



Uu(p, P + c[)6[E A (p + q) - E(p) 



00 



(2tt) 3 E(p)E A (p + q) 

where M is the number of protons or neutrons in the nucleus. Here mjv and m A are the 
masses of the struck nucleon and A, respectively, and E(p) = 



in 



N 



+ p 2 and E A (p + q) 



171 A + (p + Q) 2 are their corresponding energies; kp is the Fermi momentum and f^ v the 
inelastic single-nucleon tensor of the iV — >■ A transition. 
Introducing the standard dimensionless variables 



A 



m N 



, If 



m N 



2m 



N 



2 \2 

t = k — A , e 



2m N 



l + V 2 



and performing the angular integration in Eq. (||]) we obtain 



m N 



(4) 
(5) 



37V r* F 



SmqpK J s 



ffiv{£, Oo)ds 



where Sf — yl + rjp is the Fermi energy and 



(6) 
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e = k\ - + p 2 - \p 



T 



(7) 



the minimum energy of the struck nucleon for fixed k and A, having defined the factor 



p = i + ^(/4-i 



(8) 



which measures the inelasticity of the elementary process. 

In Eq. (H) the single- nucleon tensor /^(e, 9n) contains the angle 6 between r] and k 
given via 



cos On 



\e — rp 

KT] 



(9) 



as required by energy conservation. The condition | cos#o| — 1 then permits the response 
of the A to occur only in the range 



2k - r] F ) 2 + p\ - e F 



<A<- 



2k + i] F ) 2 + p 2 A - £ F 



(10) 



In analogy with the physics of the quasi-elastic peak [^l[| , it is convenient to introduce 
a scaling variable defined as follows 



Vl(«,A) 



£q - 1 




Ap- 1 



with £f = £f — 1- The physical meaning of the scaling variable may be deduced from 
the above equation: in terms of dimensionless variables, £f^a ^ s ^ ne minimum kinetic 
energy required to transform a nucleon inside the nucleus into a A when hit by a photon 
of energy A and momentum k. It is straightforward to check that, when = (hence 
p — 1), the ordinary quasi-elastic scaling variable ip |2l| is recovered. 

In terms of the scaling variable in Eq. (|ll|) the response region given by Eq. (|T0| ) 
simply reduces to — 1 < -0A < 1- The energy position Aap of the peak of the A response 
occurs when t/>a vanishes, namely for 



-0Z 



Aap = rp , k\ p = r(rp 2 + 1) . 



(12) 
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We now turn to a consideration of the nucleonic tensor f^, which we will obtain in 
the next section for specific N — > A currents. For any physical process this tensor must 
comply with Lorentz covariance and current conservation (Q^f^u = j t i V Q v — 0). The 
most general unpolarized second-rank tensor consistent with these requirements is |22[ 

U»(p, P + q) = -v>i(t) (g^ - + w 2 (r)V,V u - — w ^ T )E^ pa QPV° , (13) 

where w\(t), ^(t) and ws(r) are scalar functions containing the specific dynamics of the 
process and 



y, = — (p, - ^-q.) > (14) 

m N \ Q 2 J 

is a four- vector orthogonal to Q^, P^ being the struck nucleon's four-momentum. Only 
the terms involving W\ z{t) occur in EM interactions, whereas w%(t) also enters for the full 
electroweak interaction. Energy conservation via Eq. (H) implies that P ■ Q/Q 2 = — p/2: 
hence the longitudinal and transverse components of the single-nucleon tensor read 



h = foo = wi(r) + (Ap + efw 2 {r) (15) 

T 

and 



/T = /u + /22 = 2Wi(T) + 



e 2 -l 



Xe — rp " 



w 2 {t) . 



(16) 



Finally, by performing the energy integral in Eq. (||) one gets for the response functions 
in Eqs. d|D the following expressions 



Rl{k,X) 

R t (k,X) 
where 



(1 + rp> 2 (r) - w x {r) + w 2 (t)V{k, A) (1 - ^1)9(1 - ^) 



3A/~£ F /t 2 

[2^(r) + u^(r)X>(«, A)] (1 - - , 



(17) 
(18) 



Z>(«, A) 



T 
72 



(Ap + l) 2 + (Ap + l)(l + Vi)^ + 



-(l + rp 2 



(19) 
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reflects the (modest) Fermi motion of the nucleons. In fact, at the resonance peak (ip& = 0) 
Eq. ( P^D reduces to 



v(K,\) AP =tF+^ 1 e F , (20) 

and, since £f — 0.03, yields a small correction. 

2.2 Density dependence of the response functions 

In this subsection we briefly explore the density dependence of the previously deduced 
responses. This is conveniently achieved by performing an expansion in the parameter 
£p. The leading terms of the expansion of 



Rl*(k, A) = - Vi)- «t + *$&f + <kl) (21) 



are given by 



Rf = - Wl (r) + -(Ap + l)^ 2 (r) (22) 

K 

: 2 Wl (r) + {-(l + rp 2 ) + ^(Ap + l) 2 } W2 (r) , (23) 



K 

and the next terms are found to be 



= R? = - (\p + I) (I + iPI)w 2 (t) (24) 

K z 

Rf = 4 2) = ^a+v>i+v4«T). (25) 



In performing the £p —>■ limit it is of importance to realize that both responses shrink 
to the peak where k 2 = rirp 2 + 1) and A = rp. This constraint requires that T> — > when 
£f — > (see Eq. QT91)) and hence 

4 0) (k,A;^ = 0) = -w x {t) + {1 + t P 2 )w 2 {t) (26) 
4 0) (M;f F = 0) = 2w x {t) . (27) 

Expressions for the A responses in this limit will be given later. In the nucleonic sector 
one immediately obtains i?[ 0) = G\ and Rf = 2tG 2 m (cf. |2|). 

The expressions in Eqs. fll7|) , (|i~8l) and ([19|) are valid for any process involving an 
initial nucleon which is converted to an on-shell resonance of mass In particular, the 
limit p = 1 yields the response functions in the quasi-elastic peak region ]2T] . The specific 
physical process gives rise to different w\ and w 2 functions, which are evaluated in the 
next section for the iV — > A transition. 
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2.3 Nucleonic tensor 

We now evaluate the invariant functions w± and w 2 relative to the 7^ — > A process. 
These functions, being scalars, can be computed in any reference system and the most 
convenient one is found to be the rest system of the A, where the Rarita-Schwinger spinors 
take their simplest form. Let us denote with Q* = (o;*,q*) the four-momentum transfer 
in this system (the corresponding four- vector in the nucleus laboratory frame is Q^) and 
let P* = (E*,p*) be the four-momentum of the struck nucleon. The A system is then 
defined by 

pl = 0^p* = -q* = (0,0,-g*) , (28) 

having chosen the z-axis in the direction of q*. The energy conservation condition ac- 
cordingly reads 

E* = m A - u* . (29) 
Using the energy-momentum relation for the initial nucleon, namely 

m 2 N = E* 2 - p* 2 = m\- 2m A tu* + Q 2 , (30) 



we find the value of u* in terms of Q 2 

w , = ml-rn 2 N + Q 2 
2mA 

In this system has its space components parallel to q*, while its time component reads 

V °-^[ E ~ ) - ~^Q 2 ■ (32) 
Then we can easily compute the 00 and 11 components of the nucleonic tensor in 
Eq. Cp~3|) , obtaining: 



A*i = -wig n = wi (33) 

/o ° = ~ Wl [ 1 -Q 2 -) +W2 ^W 

q* 2 m 2 A q* 4 

= W ^ + W2 ^Q*- (34) 

By inverting the above equations the structure functions in the A-system are found to be 



w x = f* n (35) 

- = 5£ - P :: ) • <36) 
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Hence the problem is reduced to computing just two components, namely /q and of 
the nucleon tensor. 

With Jones Sz Scadron ]I2| we write the transition matrix element for the A excitation 
as follows 



(A\j li \N)=tfr fl0 u, (37) 
the most general form of the vertex being 

r^/j = CiT 1 ^ + c 2 r^ + c 3 r^ . (38) 

In the above the C a are (invariant) form factors and the three couplings are given byQ 

= (Q^-Qg^h 5 T 3 + (39) 

Tip = (Q p K,-Q-Kg^) l5 T 3 + (40) 

^3 



Kb = (Q?Q»-Q 2 9^HT 3 + , (41) 



where = (P^ + P A/1 )/2 and T + is the iV — > A isospin transition operator [p4 
The iV — > A tensor /* then reads 

= |— E(^i r MA«)*(4r^) , (42) 

o m N SSA 

the factor 4/3 arising from the isospin trace 

E <t\T 3 \t A ><t A \T 3 + \t>= \ (43) 
tt A 6 

and being the Rarita-Schwinger spinor describing a spin 3/2 particle. In the A rest 
frame the latter has the simple form 



u° A (0,s A ) = (44) 
«a(0,sa) = Ed^H^elwA^s') , (45) 

As' 

where the (A = — 1, 0, +1) are spherical vectors and 



u A (0,s') = 




*Note that the expressions for the transition matrix element of Dufner and Tsai and of Devenish 
et al. |Q differ from the ones of Jones and Scadron because the latter employ the set of basis vectors 
(K,Q), whereas Dufner and Devenish use (P A ,P) and (Pa,Q) respectively. Hence both ]?2 and T3 and 
the form factors C2 and C3 of these three sets of authors will be different. 
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is a zero momentum 1/2-spinor. 
Since the initial nucleon spinor is 

u(p*,s) = 



l + e* I Xs 



l +£ *'Xs 

the matrix element of the four- dimensional gamma-matrix 

Tii Ti2 

1^21 ^22 



(47) 



(4£ 



in the A-system will be related to the corresponding bispinor matrix element according 
to 

(7 ■ rf 



u A Tu 



r„ + r 



12" 



x = x' r ^x 



(49) 



'1 +e*_ 

which defines T. The general matrix element between 1/2 and 3/2 spinors can be written 

as 

v? A r u = (fsAl^rils) , (so) 

where the spin transition operator for the A 



has been introduced. 

The tensor /* evaluated in the A rest frame then turns out to be 



(51) 



[IV 



6 m A ss A 

4 rriN 

3 777a 

4 rriN 
9 UJa 



Tr 
Tr 



2r Mi r^ - ieij k T ui a k T 



(52) 



where the trace is meant to be performed in a two-dimensional non-relativistic spin space. 
By inserting Eq. ([38]) into the above expression and by performing the traces we get 



oo 



— q (E -m N ) 



9m 



N 



C 1 + ^±^C 2 + uj*Cs 2 



(53) 



8mA 



9m 



N 



(E* - m N ) C 2 (E* + m N ) 2 + C 2 (m N + m A ) 2 - C x C(m N + m A )(E* + m N ) 



(54) 
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where a form factor 



C = d + ]-C 2 {m A - m N ) + C 3 % 

2 m A + m N 



(55) 



has been defined. 

It thus appears from Eq. (^) that /q only depends upon a particular, quadratic 
combination of the three form factors Cj, namely 



G 



c 



Am A m^ 
3(m A + m N ) 



d + -(E* +m A )C 2 + u*C 3 



(56) 



commonly referred to as the Coulomb form factor. On the other hand, the expression 
in addition to the terms G\ and C 2 , also contains their cross product CC\. In the spirit 
of the familiar Sachs form factors, new form factors can however be defined in such a way 
that only squares will appear: indeed by diagonalizing the quadratic form in Eq. Q5l|) one 
obtains 



where 
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8m A , » 
9m 2 N { 



m A 



3(77i jy + TOa) 

4m n 



[G 2 M + 3G 2 E ] 



2 / ^* + mjy \ 
-m N ( C - C7i ; 

■ / 



(57) 

(58) 
(59) 

m A + m N J 

are usually referred to as electric and magnetic form factors. It is easily verified that the 
expressions in Eqs. (|56D , (|58D and ( ^9|) coincide with the definitions of the form factors 
given in [fl2 1 . Furthermore the combination G\j+3G 2 E entering in Eq. (|57D is the one which 
usually appears in the generalized Rosenbluth formula for the transverse cross section in 



-777,7V C + Ci 



?77a + ?77/v 
75* +77^ 



terms of magnetic and electric form factors p2 
Finally, by means of Eqs. (^5|,[ 



we obtain the two structure functions Wi and w 2 : 



2 

(G 2 M + 3G|) (60) 
f %G 2 C - G 2 M - 3Gl) . (61) 

\ m A / 

To finish this section it is useful to write these structure functions in an arbitrary 
system, a result which will be also used in the next section to compare with the non- 
relativistic result. Here we only consider the structure functions coming from the magnetic 



Wi 



8m a , ™ 
9m 2 N K 



m N ) 



3(m N + m A ) 
4m n 



w 2 



8Q 2 



9m A q 



*2 



(E* - m N ) 



3(m N + m A ) 



4m 



N 
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Ml nucleon multipole (proportional to G 2 M ), which is the leading contribution. Later in 
the calculations we will study the effect of including the electric and Coulomb form factors. 
In order to relate these expressions to the ones obtained using the so-called magnetic form 
of the operator, we introduce a new dimensionless form factor G defined by 

_ m N + m A 3m 2 N G M 

2m N (m N + m A ) 2 -Q 2 

Using this form factor, the magnetic invariant functions can be written 

wi = ^fq* 2 (E* + m N )G 2 (63) 
9m° N 

W ? = — ^ ' ( 64 ) 

m A q* z 

where we have used the fact that in the A-system 

(m N + m A ) 2 - Q 2 = 2m A (E* + m N ) . (65) 

In the A system we have the following expressions for the quantities involved, written in 
terms of the invariant r: 

g*2 = 4m ^ T l±^! ( 66 ) 

m N + E* = —(1 + /2 A + 2t P ). (67) 

Using these relations, the structure functions can be written in a general system as func- 
tions of T 

32 

Wl = — r(l + rp 2 )(l+/i A + 2rp)G 2 (68) 
9 



w 2 = -^ = ^r(l + /i A + 2rp)G 2 . (69) 
1 + rp z 9 



From the above equations we note that for a pure Ml transition (1 + rp 2 )w2 — W\ — 0, 
hence only the w<i structure function contributes to the longitudinal response, i.e., 

A) = ^^--w 2 (t)V(k, A)(l - Vft)0(l - ip A ) . (70) 

Finally, we end this section by obtaining the relativistic structure functions w\ and 
W2 in the limit rjp = 0. We begin with the expressions in Eqs. (|60|j6l|) written in the 
A-system which, in terms of the dimensionless variables introduced previously, read 

Mr) = ~/Wa + l) 2 r [G 2 m (t) + 3G|(r)] (71) 



A^GUr) - G 2 m (t) - 3G%(t) 



(72) 
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where = (E* — TOjv)/w7V is the nucleoli kinetic energy in the A rest frame. In the 
£f = limit the following kinematical relations hold: 



^ = ^(l + rp 2 ) (73) 

A*A 

T = — (2rp + l-/i A ). (74) 



As a consequence the single-nucleon structure functions reduce to 



w 1 (t;£ f = 0) = I( /iA + l) 2 (2rp + l- / u A )(Gi f + 3G|) (75) 



2 

u*(t;6, = 0) = + l f TP , t 1 ~ ^ -& M - 3G% ) (76) 

/ 1 -+- rp \ p A 



and the response functions in Eqs. ( p6| , p7D read 



R ( l\Zf = 0) = ^(p A + l)2(2rp + l- / i A )G 2 7 (77) 

A*A 

4°)(^ = 0) = (p A + l)2(2rp + l-/i A )(G| f + 3G|) . (78) 



3 Non-relativistic reduction of the delta current 

In this section we present a new "relativized" expression for the delta current that can be 
implemented very easily in standard non-relativistic models. Contrary to previous work 
on the iV — > A transition, where an expansion in the transferred momentum (q) and 
transferred energy (u) was performed, here we only consider expansions in powers of the 
bound nucleon momentum 7] = p/rriN, keeping the exact dependence in both k = q/2rriN 
and A = Loj2mjq. Thus, we follow closely our previous work fl"9| , [2(J where new relativized 
expressions were obtained for the electroweak nucleon and meson-exchange currents. It 
is important to realize that for high-energy conditions the traditional current operators 
- obtained assuming that also k < 1 and A 1 — are bound to fail, whereas our past 
work provides a way to include relativistic aspects into improved, effective operators for 
use with the same non-relativistic wave functions. 

Following [|rj| the full 7iVA vertex in Eq. ( ]3"%D can be decomposed into magnetic 



dipole, electric quadrupole and Coulomb (longitudinal) quadrupole contributions. In 
what follows we focus on the contribution of the magnetic term: due to the smallness 
of the E2 and C2 form factors, the magnetic term is the dominant one for both the 
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longitudinal and transverse nuclear response functions. Thus, in the particular case of 
the magnetic contribution, the ^yNA vertex current can be written in the form |T3 

r£ = ^r^(KQ) , (79) 

where we use the notation e^(KQ) = e Ptxai K a Qi with K a = (P a + P%)/2 and Q a the 
four-momentum transfer. The form factor G in Eq. (|B2J) is proportional to the mag- 
netic form factor Gm- It can be proven that the relativistic structure functions obtained 
with the current in Eq. ([79]) coincide with the general expressions given by Eqs. (60j61|) 
neglecting the electric and Coulomb form factors, i.e., Ge = Gc = 0. 

In order to proceed with the non- relativistic reduction of the ^yN — > A current, let us 
rewrite the transition matrix element for the A excitation 

(A\j^P A ,P)\N)=u A (P A ,s A )rM U (P,s) , (80) 

where u A (P A , s A ) is the relativistic Rarita-Schwinger spinor describing a spin-3/2 particle 
whose general expression is given by ||24[ 

u A (P A ,s A ) = E(f s ' 1A lf s A>e A (P A ,A) MA (P A , S ') . (81) 

As' 

Here e A (P A ,A) are the spherical vectors boosted to momentum p A and u A (P A ,s') are 
spin-1/2 Dirac spinors. Introducing the dimensionless variables rj A = ^ and e A = 
we can simply write 

e A (P A ,A) = (e A • 77 A , e A + ^^Na) (82) 



u A (P A , s') 



e A 



1 

1+£A 



Xs' ■ (83) 



Thus, following the same procedure developed in |19|, |2(| , the matrix element for a general 
four-dimensional gamma matrix T can be related to the corresponding bi-spinor matrix 
element according to 

uaTu = xl'Txs = xl> l 1 [(1 + £A)ru - a ■ r) A T 2l ] Xs + 0( V ) , (84) 

'2(1 + £ A ) 



where the expansion only to leading order in 77 has been considered. 

To make the discussion that follows easier we introduce new dimensionless variables 
which are analogous to the ones already used in the case of the nucleon quasi-elastic peak 



19|, |20|. In particular, let us define the variable A' given by the relation: e A = 1 + 2A'. It 



is straightforward to derive the following relations (valid up to first order in rj) 
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A' = — [ A 

HA 



1 

HA 



r + 



HA - 1 



( Ha - 1 
V 2 



r' + r7 • re' + 0(r? 



Ha Ha 



with /J, a as defined in Eq. ([J). The current operator T in Eq. 
to leading order in r\ in the form 



(85) 

(86) 
(87) 

can be then written up 



1 



o(v) ■ 



r = -===[(i + T')r n -<T-K'r 21 ] 

vi + T 

Moreover, in terms of the spin transition operator for the A in Eq. (]5lf) , the general 
matrix elements between 1/2 and 3/2 spinors result 



u° A Tu 



u A Tu 



(|s A |2St.KT|is) + 0(77) 



« 2 
(|S A 



I 2 

St ■ K? 

1 + t' 



2k' 



r\h) + o( v ). 



(89) 
(90) 



The non-relativistic reduction (valid up to leading order in rj) for the magnetic NA 
current operator, J M , is defined through the relation 

(ls A \J fl \±s) = {A\j fl (P A ,P)\N) (91) 

and can be obtained by using the general relations given by Eqs. (|85| - |87|) and taking into 

up 

time and space components are given by 
G 



account the fact that the magnetic current operator Tff* is diagonal in spin space. The 



Jo 
Ji 



m 



vr 



N 



St + ^2re' 
1 + t' 



e i0 (KQ) 



G 



m 



VT+r 7 { 2St • K'e oi (KQ) 



A' 



1 + t' 



eji(KQ) 



Using the relations 



e i0 (KQ) = e i0jk K j Q k ~ -2m 2 N (r) x re) 
£ji(KQ) = e jiaf3 K a QP ~ 2m 2 N e jik K k 



(92) 
(93) 



(94) 
(95) 



and taking into account the fact that re' and re are parallel vectors (see Eq. fl3?D), we 
finally get 



Jo 
J 



-2Gv / TT7 7 S t • (rj x re) + 0(rf) 
2GVT+t 1 (S ] x re) + 0(77) . 



(96) 
(97) 
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Note that the current J is of order 0(1), whereas the charge Jo is of order 0(rj). Thus, 
one expects the contribution of the NA current to be considerably more important for 
the transverse response. 

It is important to note that in order to be consistent one should treat the charge 
and current at the same level and perform an expansion of the current to order 0{rj). 



Such program can be carried out by using the techniques developed in [19], £0|. However, 
for the present case an additional simplification can be made by using the expression 
of the charge operator to order 0(r]) in Eq. (0), and taking into account the fact that, 
before performing the expansion, the original current was gauge-invariant. This invariance 
property should be valid for all the orders in the expansion (we are not expanding in q or 
uj). Hence we can relate the longitudinal component of the current with the density, i.e., 

J • k = XJq = -2GAv / TT7 7 (S t x n) ■ k . (98) 

Hence we can obtain an improved and gauge-invariant current by adding to the expression 
in Eq. (|97| ) a new piece of order 0(rj) given by 



Jgauge = -2G , Av / TT7 7 (S t X n) . 

Of course there could be additional corrections of order 0(rf) in the transverse current 
that cannot be fixed by just using the continuity equation. However, as we will show 
below, the transverse response function computed with this improved current differs from 
the exact relativistic results by terms only of order 0(t] 2 ), proving the high quality of this 
expansion for most applications in nuclear physics. Therefore the new expression of the 
current that we will consider below is 

J = 2Gv / TT^ 7 S t x (k- Xrj) . (99) 

In order to test the quality of the above expansion of the current we proceed to compute 
the response functions in a non-relativistic Fermi gas using relativistic kinematics and the 
new currents in Eqs. (p6|j99| ). We begin by computing the analytical expressions and 
compare with the exact relativistic answer. Results are shown and discussed in the next 
section. 

The calculation of the non-relativistic nucleon tensor needed to evaluate the nuclear 
response functions can be done by performing the following traces 

C = ~^[jiM=<r^ (100) 
O 777 at 

fn+f% = ^Tr[4.J T ] = 2< r , (101) 
3 m N 
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where we have introduced the longitudinal and transverse, non-relativistic structure func- 
tions and wE„ 



w 



w 



3 2>rriN 



(102) 
(103) 



In order to compare the relativistic and non-relativistic response functions it is con- 
venient to use the kinematical relations: 



1 + r' 



1 



2/i A 



(1 + fi A + 2rp) 



k 2 -2\k-t) = r(l + rp 2 ) + 0(rf 
We can then write, up to first order, the relations 



w 1 + 0(r] 2 

K 2 W\ 



+ 0( V 2 



-w 2 + 0(rf) 



r 1 + rp 2 

where w\ and w 2 are the magnetic relativistic functions given in Eqs. 
The nuclear response functions are given finally by 

37VXf 



(104) 
(105) 

(106) 
(107) 



R nr 



r>nr 
ll-r 



3A/~£f 
3Af£ F 



6(1 






V 




6(1 




T 


w 2 + 0(rf) 


V 


6(1 








6(1 




-^l)2\w 1 + 0( V 2 ) 





(108) 



(109) 



with D as given by Eq. fll^). Although we call these functions "non-relativistic", actually 
they contain enough relativistic ingredients to be high-quality approximations to the exact 
RFG result. In fact, first we use relativistic kinematics, so that the phase space and 
momentum integrals are done exactly. Second, we consider the new currents expanded to 
include effects up to order 0(rj), so that the dynamics of the problem are correct to that 
order. Comparing these expressions with the exact relativistic responses in Eqs. ( pT8| , |70D , 
we see that the relative differences between R 7 ^ and Rl, and between RVf and R? are 
of order 0(rj 2 ). In the next section we test numerically the quality of the non-relativistic 
responses of Eqs. ( |108| , |109| ). 
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4 Results 



In this section we present numerical results for the relativistic response functions in the 
region of the A-peak and check numerically the quality of our new approximation to 
the N — > A electromagnetic current. We also investigate the importance of different 
contributions in the relativistic responses. We shall present results for medium and high 
momentum transfers, ranging from q = 0.5 to 2 GeV/c, and thus will also be able to 
analyze the validity of the traditional non-relativistic calculations performed for different 
values of q. 

Before starting our analysis it is convenient to check our model with some of the 
available experimental data ||25|| . In this way we can fix some of the ingredients that enter 
in the model, in particular, the magnetic form factor Gm of the A and the modification 
of the response due to the finite A width as a consequence of its later decay into the N-tt 
channel, which we incorporate by performing a convolution with the responses for stable 
particles. 

In fig. 1 we show the transverse cross section ot for the inclusive reaction H(e, e') from 
the nucleon. The transverse nuclear cross section is defined by 

2na 2 

G T = TIC Q2~ r T , (HO) 

^ 2mjv 

where tt = j^2wi is the transverse response of a single nucleon. In this calculation 
we have included the A width by substituting for the energy-conserving delta function a 
Lorentzian shape 

5(u + m N - E A ) — ^~ (r r( f 2 ^ rM2M , (HI) 

where s = {m^ + u) 2 — q 2 is the invariant mass of the initial photon and nucleon. Results 
that include a width T(s) (which is zero at threshold and equal to the width T = 120 
MeV at resonance) are represented by solid lines in fig. 1. The dependence of T(s) on the 
invariant mass is given by [|J 

r(*) = r ^(4V , (ii2) 



where p* is the momentum of the final pion resulting from the A decay (in the A-system 
given by 

1 



(s-m 2 N -ml) 2 2 2 1 1 
~ \ — " ^nK 



113) 
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and p^ es is its value at resonance, obtained from the above expression for y/s = m&. For 
comparison, we show in fig. 1 with dashed lines results obtained by considering a constant 

width r = r . 

Other ingredients that enter in the cross section are the A form factors. We use the 
parameterization 

G M (Q 2 ) = G M (0)f(Q 2 ) (114) 
G E {Q 2 ) = G E (0)f(Q 2 ) : (115) 

where we assume that the same dependence in Q 2 is valid for the electric and magnetic 
form factors, given by the function || 

with Gg the electric form factor of the proton, for which we use the Galster parameteri- 
zation (1 + 4.97r) -2 |2(|. The above equation reflects the fact that the isobar form factor 



falls off faster than the proton form factor. Unless otherwise indicated, we take Gc = 
and use the following values [JEJ of the form factors at the origin 



G M (0) = 2.97, G E (0) = -0.03. (117) 

Later on we show the effect on the longitudinal response function introduced by consid- 
ering a C2 form factor that is different from zero. 

Our calculation (solid lines) displayed in fig. 1 is slightly below the data for the three 
values of the momentum transfer Q 2 = —0.2,-0.3,-0.4 (GeV/c) 2 , reflecting the fact 
that we have not included the background contributions of non-resonant pion production, 
which produce an additional increase of the cross section. 

In fig. 2 we show results for the nuclear inclusive cross section per nucleon from 12 C 
compared with the experimental data taken from []2], |27| . Dotted lines correspond to the 
RFG for a stable A. A more realistic model of the A peak requires the inclusion of the 
A width in the cross section, which we show with dashed lines. The nuclear responses 
including the width, Rr(q,u), are computed from the responses R(q, uj,W) for a stable 
A with mass W by a convolution 

r w max l r(W)/2 
Rv(q,u)= -7777 \" R(q,u>,W)dW , (118) 

Jm N +m, n TT (W — 171^) + T(W ) 2 / 4 

where the integration interval goes from threshold to the maximum value allowed in the 
Fermi gas model, W^ ax = (Ep + uj) 2 — (q — kp) 2 . The inclusion of the A width produces 
a broadening of the A peak and correspondingly a decrease of the strength. 
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As an illustration of how one could improve the model in the quasi-elastic-peak region, 
we also show with dot-dashed lines the quasi-elastic cross section computed with the 



PWIA model of ref. pq| . Here the mean difference with the RFG model is the inclusion 
of the momentum distribution of the finite-sized nucleus, which produces the "tails" of 
the cross section, and the binding energy of the nucleons in the nucleus, which produces a 
shift to higher energies (in the direction of data). Here for the PWIA calculation we use 
relativistic kinematics, final states are described as plane waves and the electromagnetic 
current used contains relativistic corrections to order 77. Finally, we show with solid lines 
the results computed with a hybrid model in which we add the PWIA cross section for 
the quasi-elastic contribution to the RFG result for the A contribution. 

As we can see in fig. 2, our results are below the data in the dip and A region. This 
was expected because other contributions coming mainly from two-nucleon emission and 
non-resonant pion production (not included in our model) also enter here [jl], 0, 10, 29 



However, our intention in this work is not to reproduce the experimental data nor to 
present a complete model including all of the physical contributions in this energy region, 
but instead to discuss the effect of different ingredients in the calculation and present a 
new set of improved currents specifically for excitation of the A peak that now include 
the relevant relativistic content — these could now straightforwardly be used in standard 
non-relativistic many-body models with relativistic kinematics. 

The quality of the new approximation to the relativistic, magnetic A current is shown 
in fig. 3, where the exact RFG longitudinal and transverse responses using magnetic and 
electric form factors are displayed with solid lines. Here we show just the A contribution 
to the responses. In addition we show with dashed lines the responses computed in 
the non-relativistic Fermi gas model with relativistic kinematics and the new currents 
in Eqs. ( p6| , p9D . For comparison we also show with dot-dashed lines results for the non- 
relativistic Fermi gas model and relativistic kinematics, but using the traditional non- 
relativistic current. The improvement of the description of the relativistic results using 
our currents is clear from this figure — the solid and dashed lines almost coincide. This 
proves that our expansion to order 0(rj) is precise enough to describe the A excitation in 
nuclei with negligible error for high momentum transfers. 

In fig. 4 several relativistic effects and ingredients of the calculation are analyzed. 
Therein we show the longitudinal and transverse responses for q = 0.5, 1 and 2 GeV/c. 
With solid lines we show the A peak computed within our model, while with dashed lines 
we show the A peak computed using the Peccei Lagrangian [30|. The Peccei Lagrangian 
only includes the first coupling T 1 ^ given in Eq. (^) with coupling constant C*i(0) = 2.5 
GeV -1 . This value has been chosen so that the transverse response is equal to the one 
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computed with the full Lagrangian in the A peak for q = 500 MeV/c. The importance 
of using the full vertex in Eq. (|3^) is clear from this figure. First, although the coupling 
constants can be chosen so that the transverse responses computed with both Lagrangians 
are similar for moderate q = 500 MeV/c, they begin to fail for higher g-values. These 
differences are seen to be most important in the longitudinal response, where the Peccei 
Lagrangian clearly gives an extremely large result. The reason for this unphysical behavior 
of the Peccei Lagrangian is that in Eq. (|53|) there are important cancelations among the 



C*i, C2 and C3 pieces in the longitudinal channel |TTJ. Second, as the Peccei Lagrangian is 
the one usually employed to compute the MEC contribution involving virtual A excitation 
H] , it is mandatory to use the full Lagrangian — or at least the magnetic piece — if one 
wants to compute the longitudinal contribution of the MEC in this channel. 

As reference, in fig. 4 we also show with dot-dashed lines the quasi-elastic peak re- 
sponses. For high g-values the two peaks overlap and the importance of the A in both 
responses also increases with q. This is better seen in fig. 5 where we also show the effect on 
the response functions produced by incorporating the finite A width (solid lines). Dashed 
lines correspond to results without including the isobar finite width, while dot-dashed 
lines correspond to the quasi-elastic peak. In the above results no Coulomb form factor 
has been included. In this case the contribution of the A in the longitudinal response 
is due to the Fermi motion of the nucleons inside the nucleus ||, the main contribution 
here coming from the magnetic A excitation, which is zero only for nucleons at rest. This 
is better seen in Eq. ( |100| ) where the longitudinal, magnetic single-nucleon response is 
seen to be proportional to r\\ and to the function wY which is proportional to k 2 (see 
Eq. ( |102|) ), explaining the increase with q of the longitudinal response observed in fig. 5. 
In the static limit, the longitudinal response becomes proportional to the Coulomb C2 
form factor, as shown at the end of sect. 2. 

In fig. 6 we show the dependence of the longitudinal response on the Coulomb form 
factor of the A. With solid lines we show the A peak without C2 multipoles, while the 
dashed lines include a Coulomb contribution with form factor 

G C (Q 2 ) = -0.15G M (Q 2 ) • 

We can see that the longitudinal response is quite sensitive to this form factor, especially 
for high q. This can be easily understood from the analytical expression of the structure 
functions in Eqs. ( |^D| , |5T] ) . Only the structure function w 2 depends on Gc and its depen- 
dence on the form factors is carried by the quadratic combination ^tG c — G 2 M — 3G%, 
which is not very sensitive to small values of Gc- However looking at Rl in Eq. (|T7D, 
we see that the correction due to the term w<{D is small; hence the main contribution 
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comes from the combination (1 + rp 2 )w2 — W\. The important point is that, if no C2 
term is present, that combination is exactly zero, i.e., (1 + rp 2 )w2 — w± = 0. Thus, the 
only contribution to the longitudinal response comes from the higher-order term w{D. 
This explains why the longitudinal isobar response is so small. On the other hand, if 
the C2 form factor is nonzero, there will be a contribution from the leading-order term 
(1 + rp 2 )w2 — W\ which is proportional to G 2 C . This term, although small, is of the same 
order of magnitude as the w-{D piece. Therefore, we can conclude (at least for values 
of Gc such as those assumed here) that a correct treatment of the isobar longitudinal 
response requires the inclusion of the Coulomb form factor. 

To finish this section, we have also explored the sensitivity of the response functions to 
inclusion of the electric E2 form factor, finding that both responses are quite insensitive. 
The reason is that the E2 contribution always adds incoherently to the Ml form factor in 
the combination G 2 M + 3G 2 E , as seen in Eqs. (|60| , |6T1) . Therefore, the inclusion of a small 
E2 form factor does not significantly modify any of the responses. In order to obtain 
appreciable effects due to the electric form factor, one could explore other observables 
where interferences can occur; for instance, one could analyze the angular distribution of 
the A emission in the transverse channel 



5 Conclusions 

In summary, in this paper we have obtained a new expansion of the relativistic A electro- 
excitation current to first order in T] = p/m^, maintaining the exact dependence on the 
momentum and energy transfers. We have tested this new expansion by performing a 
calculation using a non-relativistic Fermi gas model together with relativistic kinematics. 
The resulting longitudinal and transverse responses are found to be very close to the 
exact result computed with the RFG model. Therefore it is expected that the use of 
this current will provide a significant improvement when used in more realistic models of 
inclusive electron scattering from nuclei. 

We have also performed a comparison between different Lagrangians in the treatment 
of the A excitation, finding that the Peccei Lagrangian is inappropriate in the longitudinal 
channel for all of the g-values analyzed. We have studied the contribution of the isobar 
emission to the longitudinal response, in general finding a small contribution for medium 
g-values, although increasing significance as q increases. In our calculations we include 
the finite width of the A assuming a Lorentzian shape. With regards to the quadrupole 
amplitudes of the A, we have found a large sensitivity of the longitudinal response to 
inclusion of the Coulomb form factor of the isobar, especially for high q, a fact that could 
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be of importance for the present investigations of the longitudinal nuclear response. On 
the other hand, both L and T responses are found to be insensitive to the quadrupole E2 
form factor. 

Finally, it is of interest to extend the present studies to the two-particle emission 
channels which provide important contributions in this energy region. In particular, one 
expects these ideas to be relevant in the analysis of non-resonant pion production (e, e'pn) 
and two-nucleon emission (e, e'2N) reactions, where there is the hope of extracting detailed 
information on short-range correlations in nuclei. It is clear that, in order to extract 
ground-state properties from these reactions, high-g values are needed, since one wishes 
to reach a reasonably quasi-free regime in which final-state interaction effects are expected 
to be minimal. However, for such kinematical conditions, as our present studies indicate, 
relativistic effects play a role and it is important to have them theoretically under control. 
At this point in time we have explored ways to relativize the electroweak currents for 
selected situations — our intent is to continue along the same path for other related 
problems such as those mentioned above. 
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Figures 

Figure 1: Transverse cross section ap for the inclusive reaction H(e, e') in the A peak 
as a function of the A invariant mass y/s. The solid lines include a width T(s) which is 
zero at threshold. The dashed lines use a constant width T = 120 MeV. Experimental 



data are from 25 



Figure 2: Inclusive cross section per nucleon from 12 C, for two beam energies and two 
scattering angles. Dotted lines: RFG without A width; dashed: RFG including the finite 
A width; dot-dashed: PWIA for the quasi-elastic peak; solid: hybrid model obtained by 
adding the PWIA cross section in the quasi-elastic peak and RFG cross section in the 
A-peak. Experimental data are from p], |27 . 



Figure 3: Electromagnetic responses in the A-peak for 12 C with kp = 225 MeV/c, 
without a A-width. Solid: exact results within the RFG model; dashed: new expansion 
of the electromagnetic current to first order in 77; dot-dashed: traditional non-relativistic 
current. We use relativistic kinematics in all cases. 

Figure 4: Electromagnetic responses for 12 C in the RFG model with kp = 225 MeV/c, 
with A-width. Solid: A-peak Using the Jones & Scadron Lagrangian; dashed: A-peak 
using the Peccei Lagrangian; dot-dashed: quasi-elastic peak. 

Figure 5: Electromagnetic responses for 12 C in the RFG model with kp = 225 
MeV/c. Solid: A-peak including the A-width; dashed: A-peak for stable isobar; dot- 
dashed: quasi-elastic peak. 

Figure 6: Longitudinal response for 12 C. Solid: A-peak without Coulomb form factor; 
dashed: A-peak with Coulomb form factor Gc(0) = — 0.15Gm(0). 
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